PhD Chapter 2
Results
This series of files compile all analyses done during Chapter 2.
All analyses have been done with R 4.0.2.
Click on the table of contents in the left margin to assess a specific analysis.
Click on a figure to zoom it
1. Maps
1.1. General map
1.2. Parameters maps
Maps of functional traits density:
Body: non-calcified tissue
Body: calcareous
Body: calcium carbonate
Body: amorphous calcium carbonate
Body: aragonite
Body: calcite
Body: high magnesium calcite
Body: chitinous
Size: small
Size: medium
Size: large
Food: filter feeders
Food: surface deposit feeders
Food: subsurface deposit feeders
Food: grazers
Food: predators
Food: scavengers
Food: parasites
Mobility: sessile
Mobility: limited
Mobility: mobile
Lifestyle: fixed
Lifestyle: tubicolous
Lifestyle: burrower
Lifestyle: crawler
Lifestyle: swimmer
2. Rank-Frequency diagrams
We drew Rank-Frequency diagrams to study the structure of communities when considering taxa frequencies.
3. Indicators of ecosystem status
This section tests different indicators to reflect the environmental status in Baie des Sept Îles. We will consider classic methods, such as community characteristics, with functional diversity indices and other techniques. We will look at their results critically to see which could be the best for which situation.
When relevant, we used the five classes based on Environmental Quality Ratios established for the WFD and MSFD:
- 0 < status ≤ 0.2: low (red, #FF0000)
- 0.2 < status ≤ 0.4: bad (orange, #FFA500)
- 0.4 < status ≤ 0.6: moderate (yellow, #EEEE00)
- 0.6 < status ≤ 0.8: good (green, #228B22)
- 0.8 < status ≤ 1: high (blue, #0000EE)
The reference value (the denominator of the ratio) is specific to each indicator.
3.1. Specific richness
3.1.1. Methodology
We calculated a basic community characteristic, the specific richness, to see if patterns could be detected in the study area. The same calculation as for Chapter 1 have been performed for the considered stations.
Assumption: A higher species richness indicates a more desirable status.
3.1.2. Application
The following map represent the corresponding EQR status for each station:
3.2. Total density & biomass
3.2.1. Methodology
We calculated basic community characteristics, the total density and biomass of individuals, to see if patterns could be detected in the study area. The same calculations as for Chapter 1 has been performed (with the addition of biomass data) for the considered stations.
3.2.2. Application
3.3. Diversity & evenness
3.3.1. Methodology
We calculated basic community characteristics, the Shannon diversity and Pielou evenness, to see if patterns could be detected in the study area. The same calculations as for Chapter 1 has been performed for the considered stations.
Assumption: A higher diversity indicates a more desirable status.
3.3.2. Application
The following maps represent the corresponding EQR status for each station:
3.4. Taxonomic distinctness
3.4.1. Methodology
We calculated a basic community characteristic, the taxonomic distinctness, to see if patterns could be detected in the study area. The same calculations as for Chapter 1 has been performed for the considered stations.
3.4.2. Application
3.5. Functional diversity
3.5.1. Methodology
We studied functional diversity based on these species traits:
- body composition (non calcified tissue, calcareous, calcareous calcium carbonate, calcareous amorphous calcium carbonate, calcareous aragonite, calcareous calcite, calcareous high magnesium calcite, chitinous)
- body size (small, medium, large)
- food diet (filter, surface deposit, subsurface deposit, predator, scavenger, grazer, parasite)
- mobility (sessile, limited, mobile)
- lifestyle (fixed, tubicolous, burrower, crawler, swimmer)
Species were assigned to a trait using a binary code (0: absence of the trait, 1: presence). This allowed to calculate functional richness, divergence and evenness according to Laliberté & Legendre (2010).
3.5.2. Application
3.6. AZTI Marine Biotic Index (AMBI)
3.6.1. Methodology
AMBI is an ecological index that is used to detect a perturbation in an ecosystem based on its communities (Borja et al., 2000). This perturbation is linked with the organic matter loading, according to Pearson and Rosenberg (1978)’s model.
To compute this index, species are classed into five groups in relation to their tolerance to this perturbation:
- group I (GI): vulnerable species
- group II (GII): indifferent species
- group III (GIII): tolerant species
- group IV (GIV): first-order opportunistics
- group V (GV): second-order opportunistics
These groups are based on expert opinion on the physiology of species and experimental studies, but the attribution of a species to a group can be somewhat arbitrary (e.g. based on related phyla information) so it needs to be interpretated carefully. The AMBI index (also called biotic coefficient) is continuous between 0 to 6, and is calculated using this formula:
\[ AMBI = \frac{\sum_{i}^{GI-V} w_{i} . P_{i}}{100} \]
- \(P_{i}\) is the proportion of each group (percentage of the total number of species)
- \(w_{i}\) is the weighting parameter of each group (respectively 0, 1.5, 3, 4.5 and 6)
- \(i\) is the ecological group
3.6.2. Application
3.7. Multivariate AMBI (M-AMBI)
3.7.1. Methodology
M-AMBI is a complementary method that is used to calculate an Ecological Quality Ratio (EQR), a measure of the good environmental status. It is based on a multivariate ordination of the stations using the AMBI index, the species richness and the Shannon diversity. The result givises a value between 0 and 1 after comparison to reference values.
Assumption: A higher richness and diversity indicates a more desirable status.
These values are called “references” but this needs to be discussed as this vision is limited. They have been set with the 95 % percentile of the distribution. This is a recommendation by Nicolas Desroy, so that we do not detect an increase of EQR when there is a small perturbation (see work by Pearson & Rosenberg and the Intermediate Disturbance Hypothesis).
This calculation yielded 21 for S and 2.53 for H.
3.7.2. Application
The following map represent the corresponding EQR status for each station.
No clear tendancy can be discovered here, apart from the fact that the overall status seems to be “High”. Several hypothesises can explain this result:
- the M-AMBI index describes reality well, so that overall perturbation from organic matter is low
- there is a bias in the index due to the species classification in groups, originally suited for European ecosystems
- the assumptions for the reference values are not correct
- the configuration of the bay makes the perturbation small relative to the water volume and bathymetric condition
Further work is needed to determine the individual responses of somes species, along with the use of different methods to understand other perturbations and cumulative impacts.
3.8. Benthic opportunistic polychaete/amphipod ratio (BOPA)
3.8.1. Methodology
BOPA is an index that uses a relative abundance ratio of species in a community to infer a state of perturbation. Ratios with many species have been tested, and opportunistic polychaetes and amphipods have been selected to be the most pertinent (originally to detect effects of an oil-spill on soft-bottom communities, e.g. from the Sea Empress or the Amoco Cadiz). It indicates an absence of pollution when amphipods are dominant and a pollution when opportunistic polychaetes are dominant. It has been updated from its original form in 2000, and varies between 0 and \(log_{10}(2)\) (~ 0.3).
The equation is:
\[ BOPA = \left( \frac{f_{P}}{f_{A} + 1} + 1 \right) \]
- \(f_{P}\) is the relative frequency of opportunistic polychaetes (abundance / total density)
- \(f_{A}\) is the relative frequency of amphipods (abundance / total density)
We considered AMBI groups GIII to GV for polychaetes and GI for amphipods (without Jassa genera).
3.8.2. Application
These are the polychaetes and amphipods present in our species list (including the confidence score used during group classification).
| taxon_name | group | confidence_score |
|---|---|---|
| arcteobia_anticostiensis | II | 2 |
| axiothella_catenata | I | 2 |
| bipalponephtys_neotena | II | 3 |
| chone_sp | II | 2 |
| cistenides_granulata | II | 3 |
| cossura_longocirrata | IV | 3 |
| eteone_sp | III | 2 |
| euchone_sp | II | 2 |
| glycera_capitata | II | 3 |
| glycera_sp | II | 2 |
| goniada_maculata | II | 3 |
| harmothoe_sp | II | 2 |
| hediste_diversicolor | III | 3 |
| lumbrineridae_spp | II | 2 |
| maldane_sarsi | II | 3 |
| maldanidae_spp | I | 2 |
| neoleanira_tetragona | II | 3 |
| nephtyidae_spp | II | 2 |
| nephtys_caeca | II | 3 |
| nephtys_incisa | II | 3 |
| nephtys_sp | II | 2 |
| ophelia_limacina | I | 3 |
| opheliidae_spp | I | 2 |
| pholoe_longa | II | 2 |
| pholoe_sp | II | 2 |
| polynoidae_spp | II | 2 |
| praxillella_praetermissa | III | 3 |
| sabellidae_spp | I | 2 |
| scoletoma_fragilis | II | 3 |
| scoletoma_sp | II | 2 |
| scoloplos_sp | I | 2 |
| taxon_name | group | confidence_score |
|---|---|---|
| aceroides_aceroides_latipes | II | 3 |
| ameroculodes_edwardsi | I | 3 |
| ampelisca_vadorum | I | 3 |
| amphipoda | not_assigned | 0 |
| anonyx_lilljeborgi | II | 3 |
| bathymedon_longimanus | II | 3 |
| bathymedon_obtusifrons | II | 3 |
| byblis_gaimardii | I | 3 |
| caprella_septentrionalis | II | 3 |
| crassicorophium_bonellii | III | 3 |
| guernea_prinassus_nordenskioldi | III | 1 |
| hardametopa_carinata | II | 1 |
| ischyroceridae_spp | II | 2 |
| ischyrocerus_anguipes | II | 3 |
| lysianassidae_spp | I | 2 |
| maera_danae | I | 2 |
| monoculopsis_longicornis | II | 3 |
| orchomenella_minuta | II | 3 |
| phoxocephalus_holbolli | I | 3 |
| pontogeneia_inermis | II | 2 |
| pontoporeia_femorata | I | 3 |
| protomedeia_fasciata | II | 3 |
| protomedeia_grandimana | II | 3 |
| quasimelita_formosa | I | 2 |
| quasimelita_quadrispinosa | I | 3 |
The following map represent the corresponding EQR status for each station. To use EQR classification, we used the conversion method from Dauvin & Ruellet (2007).
3.9. BenthoVal index
This index is a work-in-progress by the team of Céline Labrune and Olivier Gauthier at IFREMER. This pressure score still needs to be enhanced so that more human activities are included and the score is better defined.